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Abstract 

Using an impact parameter tag to select an enriched sample of Z° — > bb events, we 
have measured the difference between the average charged multiplicity of Z° — > bb 
and Z° — > hadrons to be — flhad — 2.24 ± 0.30(stat.) ± 0.33(syst.) tracks per 
event. From this, we have derived n& — n u d s = 3.31 ± 0.41 ± 0.79. Comparing this 
measurement with those at lower center-of-mass energies, we find no evidence that 
™b — n u ds depends on energy. This result is in agreement with a precise prediction of 
perturbative QCD, and supports the notion that QCD remains asymptotically free 
down to the scale M?. 
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Heavy quark systems are a particularly good laboratory for detailed studies of 
the strong interaction and tests of the theory of Quantum Chromodynamics (QCD). 
The large quark mass Mq 3> Aqqd, where Aqcd is the QCD interaction scale, 
provides a natural cutoff in the parton shower evolution, which keeps the relevant 
space-time region compact enough to avoid the non-perturbative domain of the strong 
interaction. Recently it has been recognized that, within the context of perturbative 
QCD, this cutoff allows a stringent constraint to be placed on the difference in light 
hadron production between e + e~ annihilation into heavy and light quarks [1]. In 
particular, it is expected that to 0([a s (W 2 )] 1/2 (MQ/W 2 )) (~ 0.1 track at W = M z ), 
the difference between the total mean charged multiplicity in light quark (q = u, d, s) 
events and the mean charged multiplicity of radiated 'non-leading' hadrons in heavy 
quark (Q = b, c) events, excluding the decay products of the 'leading' long-lived 
heavy hadrons, should be independent of center-of-mass (cms) energy W. This is 
a striking prediction, in that the total multiplicity is known to grow faster than 
logarithmically with W. Furthermore, to 0(a s (MQ)n u d s (MQ)) (~ 1.2 tracks for 
Q = b), this multiplicity difference should be equal to n u d s {\feMQ) , the mean charged 
multiplicity for e + e~ annihilation to light quarks at the reduced cms energy ^/eMg, 
where In e = 1. A test of this hypothesis provides the opportunity to verify an accurate 
prediction of perturbative QCD, and to probe the validity of perturbative calculations 
down to the scale Mq. In addition, this hypothesis is in direct contradiction with 
the hypothesis of flavor-independent fragmentation [2,3], which suggests that the 
non-leading multiplicity associated with heavy quark production at a given cms energy 
W should be equal to the total light quark (u, d, s) event multiplicity at the reduced 
cms energy (1 — (xq))W, where xq = 2 • Eq/W is the heavy hadron energy fraction 
after fragmentation. 
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Recent tests of these hypotheses [1,4] made use of a measurement of the mean 
charged multiplicity of Z° — ■> bb events from the statistically-limited data sample of 
the 1990 run of the Mark II detector at the SLAC Linear Collider (SLC), and were not 
able to demonstrate a clear preference for either. Here, we present a more accurate 
measurement based on the 1992 run of the SLC Large Detector (SLD) experiment, 
during which a total of 420 nb^ 1 of electron-positron annihilation data were recorded 
at a mean cms energy of 91.55 GeV. 

The SLD is a multi-purpose particle detector and is described elsewhere [5]. 
Charged particles are tracked and momentum analyzed in the Central Drift Chamber 
(CDC), which consists of 80 layers of axial or stereo sense wires in a uniform axial 
magnetic field of 0.6T. In addition, a silicon vertex detector (VXD) [6] provides an 
accurate measure of particle trajectories close to the beam axis. With the exception 
of the hadronic event trigger, this analysis relied exclusively upon the information 
from these two tracking systems. 

While the multiplicity measurement relied primarily on information from the 
CDC, the more accurate impact parameter measurement provided by the addition 
of the VXD information to the CDC tracks was used to select a sample enriched in 
Z° — > bb events. All impact parameters used in this analysis were for tracks projected 
into the plane perpendicular to the beam axis, and were measured with respect to 
an average primary vertex (PV) derived from fits to events close in time to the event 
under study. The impact parameter d was derived by applying a sign to the distance 
of closest approach such that d is positive when the vector from the PV to the point 
at which the track intersects the thrust axis [7] makes an acute angle with respect 
to the track direction. Including the uncertainty on the average PV, the measured 
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impact parameter uncertainty for the overall tracking system approaches 15 jim 
for high momentum tracks, and is 80 fim at p^VshiO = 1 GeV/c, where pj_ is the 
momentum transverse to the beam axis, and 6 the angle relative to the beam axis. 

Events were classified as hadronic decays of the Z° provided that they contained 
at least 7 tracks which intersected a cylinder of radius tq = 5 cm and half-length 
zq = 10 cm surrounding the average PV, a visible charged energy of least 20 GeV, 
and a thrust axis satisfying | cos 9thrust\ < 0.7. The resulting sample contained 5449 
events. Backgrounds in this sample were estimated to be ~ 0.1%. 

For the purpose of multiplicity counting, a loose set of requirements was placed 
on reconstructed tracks, while stricter requirements were placed on tracks used 
to measure impact parameters. 'Multiplicity quality' tracks were required to: i) 
have j9j_ > 0.12 GeV/c; ii) have | cos 6*| < 0.8; and iii) intersect a cylinder of 
(ro, zq) = (1.5, 5.0) cm. 'Impact parameter quality' tracks were required to: i) have 
| cos ^| < 0.8; ii) intersect a cylinder of (ro,zo) = (0.3,1.5) cm; iii) have at least 
one VXD hit; iv) have ad < 250/im; and v) have x 2 /d.o.f. for the CDC-only and 
combined CDC/VXD fits of less than 5.0 and 10.0, respectively. 

A Z° — > bb enriched sample was selected by dividing each event into two 
hemispheres separated by the plane perpendicular to the thrust axis, and requiring 
two or more impact parameter quality tracks in one hemisphere with normalized 
impact parameter d/ad > 3.0 [8]. Restricting the tag to tracks from a single 
hemisphere allowed potential tagging bias to be reduced by measuring the multiplicity 
in the hemisphere opposite to the tag. Monte Carlo (MC) studies indicate that 
this tag is 50% efficient at identifying hemispheres containing B hadrons in selected 
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hadronic events, while providing an enriched sample of 72% purity. The tag selected 
1829 hemispheres. 

In determining the total charged Z° — > bb multiplicity n b , we minimized 
systematic errors by measuring 5n b = n b ~ nhad, and then adding back in the total 
hadronic charged multiplicity n^, which has been accurately determined by other 
experiments [9]. In terms of the uncorrected mean reconstructed multiplicities 7% 
(m t ) of the total hadronic (hemisphere opposite tag) samples [4], 

5n b = (1 - R b ){n dk + n n i - n udsc ), 

where n n i and n u dsc satisfy 

rTTh = Ch,udsc{^- - Ph)n u dsc + C^dkPh^dk + Ch,nlPhn~nl, 

2mt = C t ,udsc(l - Pt)n u dsc + C^dkP^dk + C t ^iP t n n h 

and where and Pt are the fraction of Z° — > bb events in the hadronic and tagged 
samples, determined by MC studies to be 0.223 and 0.724, respectively. We have 
used the Standard Model value R b = T{Z° -> bb)/T(Z° -> hadrons) = 0.217 [10]. 
We have separated the Z° — > bb multiplicity into two components, one associated 
with the decay of the B hadrons (dk), and one associated with the remaining 
non- leading system (nl), in order to take advantage of measurements from the T45 
which constrain both the multiplicity and spectrum of B hadron decay products 
[11,12]. Here = 10.88 ± 0.22 is twice the B hadron decay multiplicity from 
the T45 [11], with an additional uncertainty of ±0.10 tracks included to account for 
the uncertainty in the production fractions and decay multiplicities of the B s and 
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B baryons. The constants Cij account for the effects of detector acceptance and 
inefficiencies, and biases introduced by the event and tagged sample selection criteria. 
The Cij were evaluated, using a MC simulation of the detector, as the ratio of the 
number of multiplicity quality tracks to generated charged multiplicity tracks for the 
six sub-samples. We have included in the generated multiplicity any charged track 
which is prompt, or is the decay product of a particle with mean lifetime less than 
3 x KT 10 s. 

Because of the exclusion of tracks with very low momentum or large | cos 6* | , 
the constants Cij are somewhat dependent on the model used to generate MC 
events; we have used JETSET 6.3 [13] with parameter values tuned to hadronic e+e~ 
annihilation data [14]. The resulting values for the Cij were 0.855, 0.905, and 0.810 
for C hjUdsc , C M fc, and C hjnh and 0.870, 0.904 and 0.818 for C t , u dsc, C t ,dk, and C t , n h 
respectively. 

The uncorrected mean charged multiplicity for all hadronic events was found to 
be = 17.29 ± 0.07 tracks, while the mean charged multiplicity opposite tagged 
hemispheres was found to be m t = 9.28 ± 0.09 tracks. Combining these values with 
the Cij via the above relations yields Snj, = 1.94 ± 0.30(stat.) tracks. 

We have investigated a number of systematic effects which may bias the measured 
value of 5ni,. Dividing by the overall reconstruction constant Ch )U dscb — 0.855 
provides a measurement of the total hadronic multiplicity Jihad — 20.21 ± 0.08(stat.). 
This value is lower than the world average 20.95±0.20 [9], indicating that the detector 
simulation overestimates the mean SLD tracking efficiency by ~ 3.5%. We account 
for this by reducing all reconstruction constants Cij by this amount, leading to a 
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correction of +0.10 ± 0.10 tracks in 5n\y We have conservatively set the systematic 
error in the correction to be equal to the size of the correction itself. 

After correcting for overall tracking efficiency a comparison of the p± distribution 
between data and MC shows good agreement for the untagged sample, but an excess 
of ~ 15% for data tracks opposite tagged hemispheres with pj_ between 0.12 and 0.50 
GeV/c, accounting for ~ 3% of all reconstructed tracks in this sample. Since there 
are currently no empirical constraints on the p± distribution of non-leading tracks 
in Z° — > bb events, we have assumed that this excess is due to improper modelling 
of the non-leading tracks by the JETSET MC, which to this point has been tuned 
only to the global features of inclusive Z° — > hadrons data. We compensate for this 
discrepancy by applying a further correction to dn^ of +0.20 ± 0.20 tracks, where 
again we conservatively assign an uncertainty equal in magnitude to the correction. 
In addition, we have studied the behavior of Snj, when numerous other experimental 
parameters, such as tracking and event selection requirements, were varied over wide 
ranges. As a result of these studies, we assign an additional systematic uncertainty 
of ±0.15 tracks due to the uncertainty in charged-particle spectra modelling. 

We have compared the fraction of tagged hemispheres ff ata = 1829/10898 = 
0.168 + 0.004 to the MC expectation jf AC = 0.157, assuming the world average value 
of Rb = 0.220 ± 0.003 [15]. If we conservatively assume that this difference is due 
entirely to extra Z° — > udsc contamination in the tagged sample, the corresponding 
change in Srtb is 0.21 tracks. Since impact parameter reconstruction errors tend to 
produce correlated changes in the Z — > udsc and Z° — > bb tagging efficiencies, the 
true uncertainty is somewhat less than this. From MC studies of tracking errors 
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which produce the observed difference in ft, we estimate the systematic error due to 
the tagged sample purity to be ±0.15 tracks. 

An additional systematic error of ±0.12 tracks arises from limited MC statistics. 
Combining these uncertainties in quadrature, and including the two corrections 
discussed above, we find 

5n b = 2.24 ± 0.30(stat.) ± 0.33(syst.) tracks. 

The effects of initial state radiation, and the ~ 0.2 GeV difference between the mean 
cms energy of 91.55 GeV and the Z° peak, are small, and no correction has been 
made. Adding back in the world-average total hadronic multiplicity at the Z° peak 
n had = 20.95 ± 0.20 [9] then yields 

n b = 23.19 ± 0.30(stat.) ± 0.37(syst.) tracks. 

To test the energy independence of the difference between the total multiplicity 
in light quark events and the non-leading multiplicity in Z° — > bb events, we make 
use of lower cms energy measurements of the e + e~ — > bb multiplicity from the 
PEP and PETRA storage rings. Assuming the energy independence of the decay 
multiplicity of B hadrons produced in e + e~ annihilation, it is equivalent to test the 
quantity An b = n b — n u d s - Results for this quantity for the various lower cms energy 
experiments are summarized in Ref. [1]. Applying the procedure presented in Ref. [1] 
to the SLD measurement to remove the contribution from Z° — > cc, we arrive at the 
result 

An b = 3.31 ± 0.41(stat.) ± 0.53(syst.) ± 0.58(n c ) tracks, 
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where we have constrained n c to lie between n u ^ s and 715, yielding n c = 21.9 ± 2.0 
tracks. 

Figure 1 shows rihad an d An^ as functions of cms energy. The Artf, data, with 
the additional lever arm provided by the SLD measurement, are seen to be consistent 
with the hypothesis of energy independence, in marked contrast to the steeply rising 
total multiplicity data [16]. Due to differing measurement techniques, results for An^ 
at PEP/PETRA energies are largely uncorrelated with those at the Z° peak. A 
linear fit to the An& data yields a slope of — 1.0 ± 1.1 tracks/In (GeV), consistent with 
at 0.9 standard deviations. Also shown is the perturbative QCD expectation for 
the value of An&. Averaging the SLD result with previous measurements [1], we find 
that An^ omb = 3.83 ± 0.63, within 1.1 standard deviations of the perturbative QCD 
expectation of 5.5 ± 0.8 ± 1.2(theory) [1]. 

The hypothesis of flavor-independent fragmentation [2,3], which provides that 
nb(W) — fldkiW) = n u rf s ([l — (xq)]W), implies that An;, decreases with cms energy in 
proportion to n u ds(W) [1], in contradiction with the perturbative QCD expectation. 
Figure 2 shows a comparison between non-leading multiplicity n^iW) —ndk(W) and 
?Ws([l — ( x q)]W), as a function of non-leading energy [1 — (xq)]W. When the 
SLD result is included, a linear fit to the residuals (Fig. 2b) yields a slope of 
s = 1.91 ± 0.65 tracks/ln(GeV), inconsistent with the hypothesis of identical energy 
dependence (s = 0.0) at the level of 2.9 standard deviations. 

In conclusion, we have measured the difference in the mean charged multiplicity 
between Z° — > bb and Z° — > hadrons to be <5n& = 2.24±0.30(stat.)±0.33(syst.) tracks 
per event, from which we calculate the multiplicity difference between Z° — > bb and 
Z° -> uds to be An b = 3.31 ± 0.41 (stat.) ± 0.53(syst.) ± 0.58(n c ) tracks. Comparing 
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our measurement with similar results from lower energy e + e~ annihilation data, we 
find no evidence that An;, depends on cms energy. This energy independence is 
in agreement with the precise perturbative QCD expectation, and indicates that 
QCD remains asymptotically free down to the scale Mjj. Our measured value 
is in reasonable agreement with the less precise QCD prediction that An& = 
n<lk — nuds(VeMq). Including our measurement, the cms energy dependence of the 
non-leading multiplicity in e + e~ annihilation to b quarks is inconsistent with that 
of the hypothesis of flavor-independent fragmentation at the level of 2.9 standard 
deviations. 

We thank the personnel of the SLAC accelerator department and the technical 
staffs of our collaborating institutions for their outstanding efforts on our behalf. We 
also thank Valery Khoze for helpful and motivating discussions. 
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Figure Captions 



Figure 1. Energy dependence of the total multiplicity [16] and the multiplicity 
difference An;, [1,16] between e + e _ — > bb and e+e~ — > -urfs events. A linear fit 
to the energy dependence of An^ yields a slope of s = — 1.0 ± 1.1, consistent with 
the hypothesis of energy independence (s = 0.0). The horizontal lines are the 
expected value and la range for Arl^ = n^fc — n u ds(v^^&)> given by lower-energy 
total multiplicity data in accordance with perturbative QCD (see text). 

Figure 2. a) Non-leading multiplicity n n i = — n^k in e + e~ — > 66 vs. non-leading 
energy (1 — (x^))lF [16]. The solid line is a fit [4] to e+e~ — > -urfs multiplicity as a 
function of W. The error on this fit (dotted lines) is dominated by the uncertainty on 
the removal of the heavy quark {Q = c, b) contribution to the measured nhad(W). b) 
Residuals of a). 
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